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Abstract 

A class of simple superposition models based on the Glauber picture of multiple 
collisions is compared with the data on the centrality dependence of the multiplicity 
distributions in a central rapidity bin. We show how the results depend on the 
specific assumptions concerning the distributions in the number of participants and 
their relations to the distributions of the number of produced hadrons in various 
phase space bins. None of the versions of the model describes satisfactorily the 
centrality dependence of the scaled dispersion. 
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1 Introduction 



In a recent note [T] we have analyzed the PHENIX data j2] on the centrality dependence 
of the multiplicity distributions in a central pseudorapidity bin for heavy ion collisions. 
We have shown that a superposition model which presents the final state as a simple 
superposition of the final states from nucleon-nucleon collisions (described by the PYTHIA 
8 generator [3]) fails to describe the data. Even if the average number of nucleon-nucleon 
collisions is treated as a free parameter to be fitted to the average multiplicity for each 
centrality bin, the dispersion values for these bins disagree with the data. 

We have suggested two possible reasons for this failure. First, the multiplicity distribu- 
tions for the pp collisions measured by PHENIX are not reproduced correctly by PYTHIA 
with the default parameter values; the fluctuations in data are significantly bigger than in 
the model. Second, it is well known that the average multiplicities in heavy ion collisions 
are not exactly proportional to the number of "wounded" nucleons |1] in colliding nuclei 
[5J . Thus it does not seem reasonable to assume that the final state may be considered as 
a simple superposition of final states from nucleon-nucleon collisions. 

One may try to remove the deficiencies listed above taking into account two following 
remarks. First, it is well known that the multiplicity distributions in any pseudorapidity 
bin for the pp collisions may be successfully parametrized by negative binomial distribution 
[5] . Such a distribution with the parameters fitted to the data should be thus used instead 
of that from PYTHIA with the default parameter values. Second, the wounded nucleon 
model is known to work much better when modified to the "wounded quark model" [7]. In 
a simple version of this model presented recently [8] each nucleon consists of a quark and 
a diquark. The nucleon-nucleon collision may be well approximated by the interaction of 
just one component from each nucleon (both contributing similarly to the multiplicity of 
the final state). However, for the nucleon interacting more than once in the heavy ion 
collision, both components are likely to interact and one gets a double contribution to the 
multiplicity of the final state. 

Another possible improvement of the model is the inclusion of fluctuations in the 
number of participants for a given impact parameter. The reliable description of such 
fluctuations requires the detailed knowledge of the nuclear structure, but one may estimate 
the influence of this effect by considering simple distributions "bracketing" the physically 
plausible distributions from below and above. 

However, such modifications are difficult to be built into the full-fledged Monte Carlo 
generator (e.g., PYTHIA). Therefore we construct a class of simple generators which, for 
consecutive values of the impact parameter, produce the samples of "events" consisting 
only of two numbers. These are the numbers of charged hadrons in the phase space regions 
corresponding to the two detectors used in the PHENIX experiment. Such "events" are 
then processed exactly in the same way as the PHENIX data. 

In the next section we describe the construction of such generators, introducing grad- 
ually the improvements discussed above and comparing the consecutive versions with the 
data. In the last section the results are summarized and some conclusions are drawn. 

2 Simple models and the PHENIX data 

The PHENIX data for the multiplicity distributions in a small central pseudorapidity bin 
are collected for various beams, energies and the centrality classes defined by multiplicity 
ranges in the auxiliary "BBC counters" covering the pseudorapidities r\ from the range 
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3 < \r)\ < 3.9. As before [lj, we concentrate here on the AuAu data for 200 GeV, as the 
dependence on the energy and atomic number is not very difficult to reproduce. 

To implement the superposition idea we use the code based on the Glauber formalism 
[9] for the heavy ion collision which calculates for each impact parameter the number of 
interacting nucleons in two colliding nuclei [TO]. This code was already used for comparison 
of the wounded nucleon and the wounded quark models [11]. In the following we assume 
that the average multiplicities in the BBC counters < ubbc > and in the central counter 
< n c > are proportional to the global number of interacting nucleons in both nuclei N p : 

< n BBC >= aN p , < n c >= f3N p . 
2.1 Wounded nucleons 

Let us assume a simple geometric distribution of the impact parameter, in which the 
number of events is proportional to b. To relate the generator to the experimental data we 
need the values of the two constants a and 0. To reproduce the experimental distribution 
of Ubbc presented by PHENIX [12] we take a = 5.2. A fit to the dependence of < n c > 
on N p shown in ref. [2] gives the value of /3 = 0.18. 

In the PHENIX data each centrality class (defined by the limits on Ubbc) was labeled 
by the value of N p corresponding to this range of ubbc i n the HIJING event generator. 
In our model this correspondence is only slightly different; the average values of N p in 
all classes (defined by the same bounds on tlbbc) here and in the next subsections are 
shifted in the worst case by a few percent. For consistency, in what follows, we use our 
values of N p on the x axes of the plots. 

The simplest assumption on the multiplicity distributions in the BBC and central 
counters (for fixed b and N p ) is to describe them by the Poisson distribution. Then one 
may calculate the scaled dispersion of the distribution in n c (for the central counter) for 
each centrality class as defined by a cut in Ubbc- More precisely, we calculate for each 
value of b the number of the events N ev to be generated (proportional to b) and the values 
of N p , < ubbc > and < n c >. Then we generate N ev values of ubbc according to the 
Poisson distribution with average < tlbbc > and register the numbers N l ev of the values 
falling in the ranges which define consecutive centrality classes^]. Finally, these values 
of N l ev are used to decide how many values of n c should be generated according to the 
Poisson distribution for each class. The final distribution of n c for each class is the sum 
of the distributions generated for all values of b. 

This procedure by definition reproduces correctly the experimental dependence of the 
value of < n c > on the corresponding average value of N p defined for the consecutive 
centrality classes. This is illustrated in Fig.l, where for transparency the big experimental 
error bars are omitted. The real test of the model is thus the centrality dependence of 
the scaled variance 

lu = (< n 2 c > - < n c > 2 )/ < n c > . 

The results are shown and compared with the PHENIX data in Fig. 2 together with the 
results from two other versions of the model (to be discussed later). Here and in the 
following the errors are comparable with the size of the data points. 

3 An independent consistency check of our procedure is the fact that the numbers of events in all classes 
defined by the limiting values of ns bc are almost the same (with a few percent accuracy) ; remember that 
these values were selected by PHENIX to divide the global sample of events into the equally populated 
bins. We repeat this check in the following for all versions of the model. 
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Fi gure 1: The average multiplicity in the central detector for the PHENIX pp and AuAu data (+ marks) 
and the superposition model (a; marks) . The results of other versions of the model, shown as squares and 
stars, will be commented upon later 



We see that the fluctuations are strongly underestimated and increase monotonically 
with N p in disagreement with the data. In fact, the values are quite similar (although 
slightly lower) to the results obtained by superimposing the pp events generated from 
PYTHIA [lj. This is so because the Poisson distribution is only slightly narrower than 
the distribution predicted by PYTHIA for the pp collisions. 

In the next step we replace the Poisson distribution by the negative binomial distri- 
bution (NBD): 
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with parameters fitted to the PHENIX pp data. More precisely, if we assume that a 
single participant nucleon yields (for some part of the phase space) a charged hadron 
distribution described by the NBD with the parameters < n > and k, we expect for the 
pp collisions the NBD with parameters 2 < n >, 2k, and for a heavy ion collision with 
N p participants the NBD with parameters N p < n > and N p k. The rest of the procedure 
remains unchanged. 

There is an uncertainty connected with this prescription for the NBD parameters. 
The PHENIX pp data for the central bin are collected on condition that there is at least 
one charged hadron in each of the two BBC counters. Thus we cannot assume that the 
exact values of the measured parameters of this distribution should be used to predict the 
distribution for given value of N p , where obviously not all the participants must produce 
the hadrons falling into these counters. Therefore in the following we compare the results 
for different values of the parameters for the "elementary" NBD in n c (and Ubbc)- 

We have performed the calculations with the average multiplicities of the n c and ubbc 
distributions calculated for each N p as before and with three different choices of the k 
parameter. In the first calculation, for all values of N p we have taken approximately the 
same ratio of k c to < n c > as in the PHENIX pp data (about 6) and assumed that for the 
distribution of ubbc this ratio is the same. In the second one, we assumed much broader 
distribution of ubbc, with kj < n >~ 1. Finally, we have chosen the distribution of n c 
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Figure 2: The scaled dispersion for the PHENIX pp and AuAu data (+ marks) and the superposition 
model with Poisson distributions (x marks) and NBD distributions with three choices of the k parameters 
(stars, open and full squares) 



with the ratio of k c to < n c > smaller than in the PHENIX pp data by a factor of 0.5. 
The last two choices may be regarded as maximizing the fluctuations due to the spread 
of multiplicities which define the centrality classes. In fact, the third choice gives such 
a broad distribution of n c for a given N p that it may be regarded as breaking the basic 
assumption of our class of models: building the heavy ion collisions from the elementary 
pp collisions. Indeed, k c j < n c >~ 3 corresponds to the value of to — 1 twice as big as in 
the PHENIX pp data. 

The first two choices give very similar results, significantly higher than the Poisson dis- 
tributions, but the fluctuations in heavy ion collisions are still underestimated (see Fig.2). 
The scaled dispersion increases monotonically with the number of participants, contrary 
to the data. For the third choice we match the data for a low number of participants, 
but not for the central events, where the model results keep increasing while the data fall 
down. 



2.2 Wounded quarks 



In this subsection a more refined model based on the "wounded quark" idea [7], [8j is 
developed. We shall modify the code [10] which was used to calculate for each value of 
the impact parameter b the global number of participants in two colliding nuclei 
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Similar arguments allow to calculate the global number of nucleons which interacted 
exactly once 
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and the number of those, which interacted at least twice 
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with the obvious condition N p = N p + N p . Assuming that the multiple interaction results 
in doubling the average multiplicity and the value of the k parameter one gets for each 
value of b the multiplicity distribution given by 
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where 7 = Np/(Np + N%). Note that the average multiplicity in this distribution is 

< n >= [jN* + 2(1 - 7 )iV p 2 ] <n >. 

To reproduce the dependence of < n BB c > and < n c > on < N p > we needed the values 
of < no > of 2.7 and 0.095, correspondingly. The ratio k/ < no > was assumed about 6, 
as before. Now we generate the distributions in the BBC and central counters according 
to this distribution and proceed as before. The same procedure is also performed for the 
Poisson distribution. The dependence of the average multiplicity on N p , shown in Fig.l as 
squares, is exactly the same as in the previous version of the model (shown as the x-signs). 
The results for the scaled dispersion are shown in Fig. 3. The values are slightly lower, 
but the main feature of Fig. 2 remains unchanged: the dispersion increases monotonically 
with the number of participants, in disagreement with the data. 
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Figure 3: The scaled dispersion for the PHENIX pp and AuAu data (+ marks) and the "wounded 
quark" model with the Poisson distributions (x marks) and NBD distributions (stars) 



Thus, contrary to our speculations formulated in [lj, the replacement of the "wounded 
nucleon" by the "wounded quark" model did not improve the situation. The reason is 
that the fluctuations of n c for a given centrality class result mainly from the fluctuations 
in the number of sources contributing to the values of tibbc i n the range defining this 
class. It does not really matter if these sources are nucleons or quarks (diquarks); for 
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the linear relation between their number and average value of Ubbc the results are quite 
similar. The slight reduction of the dispersion values is probably a consequence of the 
fact that the number of sources is now bigger. 

The values of the scaled dispersion may be increased by decreasing the k c / < n c > 
value, as in the previous subsection. However, the centrality dependence remains mono- 
tonic in disagreement with the data. 



2.3 Fluctuating number of participants 

There is still one assumption in the simple model we use which may be modified without 
abandoning the basic superposition idea. Using the code pi)] one can calculate a well 
defined value of N p (or and N£) for each value of b. Obviously, in the real collision 
of nuclei this number may fluctuate and the calculated value should be treated rather 
as an average. Moreover, there is a simple reason why such fluctuations should not 
increase monotonically with centrality: for most central collisions almost all nucleons 
are "wounded" and the average is close to the maximal possible value. This may damp 
the fluctuations and eventually may lead to non-monotonic dependence of the scaled 
dispersion on the number of participants. 

In what follows we consider different distributions which may be assumed for N p 
around this average. The reliable choice would require a thorough knowledge of the 
nuclear structure, and in particular of the two nucleon correlations in the position space. 
However, one may use simple lower and upper limits to the fluctuations. As the first 
one, we use a simple binomial distribution of the number of participants with the average 
< N p > calculated from the code [10]; its maximal value is obviously given by 2 A. Let us 
remind here that the dispersion squared of this distribution is given by 

D 2 =< N p > (1- < N p > f2A). 
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Figure 4: The scaled dispersion for the PHENIX pp and AuAu data (+ marks) and the superposition 
model with the binomial smearing in the number of participants and the Poisson distributions (x marks) 
and the NBD distributions (stars) 



Now the procedure is more involved. For each value of the impact parameter b we 
calculate < N p > as before, then generate randomly the values of N p from the binomial 
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distribution with this average, calculate the average values of tibbc an d n c (proportional 
to N p ) and generate the values of these variables according to the Poisson or NBD distri- 
butions. The results are shown in Fig. 4. 

Contrary to the naive expectations, the fluctuations in the number of participants 
did not increase the values of the scaled dispersion. In fact, the values for the NBD 
distributions are significantly lower than in Fig.2 for the same values of the k parameter. 
This may be interpreted as the result of the weaker correlation between the range oIubbc 
and of n c , since these two values are generated independently and the range of N p is now 
broader for each "class of centrality". 

The second choice, maximizing the spread of N p for a given b, is just a flat distribution. 
More precisely, we assume that for the < N p > value below Aj 4 the probability of any 
value of Np between and 2 < N p > is the same, for the < N p > value between Aj '4 and 
3A/A the allowed range is between < N p > -A/A and < N p > +A/A and for < N p > 
above 3 Aj 'A the allowed range of N p is between 2A — 2<N P > and 2 A. Now the dispersion 
squared is given in these three ranges of < N p > by 

2 < N p > (2 < N p > +l)/6, A(A/2 + l)/2A, (2 A- < N p >){AA - 2 < N p > +l)/6. 
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Figure 5: The scaled dispersion for the PHENIX pp and AuAu data (+ marks) and the superposition 
model with the fiat smearing in the number of participants and the Poisson distributions (x marks) and 
the NBD distributions (stars and squares) 



These values are obviously much higher than in the case of the binomial distribution 
considered before. Repeating the procedure described above we get the scaled dispersion 
as shown in Fig. 5. The results for the average multiplicity are shown as stars in Fig.l. 

For the Poissonian distribution and for the NBD with realistic values of kj < n > the 
results are almost the same as in Fig.4. If we choose the value of kj < n > much lower 
than in the pp data, as already used in Fig.2, we get non-monotonical dependence of u> 
on < N p >. However, the shape disagrees with the data and the values are still much too 
low. Thus we conclude that by increasing the fluctuations in the number of participants 
we are not able to reproduce the data. 
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3 Conclusions and outlook 



We considered the multiplicity distributions in the central rapidity bin in the heavy ion 
collisions for various "centrality classes" defined by the multiplicity in another rapidity 
bin. A class of simple generators is constructed basing on the assumption that the final 
state is a superposition of states obtained separately from each participant nucleon. We 
considered the modifications resulting from counting the "wounded quarks" instead of the 
"wounded nucleons" and from the fluctuations in the number of participant nucleons for 
the given value of the impact parameter. None of the versions of the superposition model 
considered is compatible with the data. In particular, the observed slow increase of the 
scaled dispersion at moderate centralities and the decrease for the most central events is 
not reproduced. 

The reason why the modifications introduced did not yield the expected improvement 
seems to be the particular definition of "centrality classes" used in the PHENIX experi- 
ment. Both the "wounded quark" idea and the fluctuations in the number of participants 
may give saturation (or even decrease) of the fluctuations for highest centrality assuming 
that the centrality is defined by a fixed range of the number of participants. However, 
when it is defined by the multiplicity of hadrons in some detector (as in the PHENIX 
data) these expectations are not justified, as we have seen in the model calculations. 

Our conclusions are that it seems to be difficult (if not impossible) to describe the 
dependence of the scaled dispersion on centrality measured by PHENIX in the framework 
of superposition models. This suggests the presence of some collective effects in these 
data. 
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